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Abstract. This paper presents an overview of the radiative transfer problem of calculating the 
spectral line intensity and polarization that emerges from a (generally magnetized) astrophysical 
plasma composed of atoms and molecules whose excitation state is significantly influenced by 
radiative transitions produced by an anisotropic radiation field. The numerical solution of this non- 
LTE problem of the 2nd kind is facilitating the physical understanding of the second solar spectrum 
and the exploration of the complex magnetism of the extended solar atmosphere, but much more 
could be learned if high-sensitivity polarimeters were developed also for the present generation 
of night-time telescopes. Interestingly, I find that the population ratio between the levels of some 
resonance line transitions can be efficiently modulated by the inclination of a weak magnetic field 
when the anisotropy of the incident radiation is significant, something that could provide a new 
diagnostic tool in astrophysics. 
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1. INTRODUCTION 

The non-LTE problem of the 1st kind is that described in Mihalas's (1978) book [Q]]: it 
consists in calculating, at each spatial point of the astrophysical plasma model under con- 
sideration, the values of the atomic level populations that are consistent with the intensity 
of the radiation field generated within the medium. This requires solving jointly the ra- 
diative transfer (RT) equations for the specific intensities and the statistical equilibrium 
equations for the level populations. Nowadays, this non-local and non-linear problem is 
routinely solved via the application of "clever" iterative schemes, where everything goes 
as simply as in the A-iteration method, but for which the convergence rate is extremely 
high [e.g., the reviews by^H]. The success of this international meeting to honor Prof. 
Dimitri Mihalas on the occasion of his 70th birthday is a clear demonstration that the 
numerical solution of this multilevel radiative transfer problem has enabled many im- 
portant advances to be made in 20th century astrophysics, including the modeling of the 
Fraunhofer spectrum and the development of the field of astrophysical spectroscopy. 



Invited review paper published in Recent Directions in Astrophysical Quantitative Spectroscopy and 
Radiation Hydrodynamics. Proceedings of the International Conference in Honor of Dimitri Mihalas for 
His Lifetime Scientific Contributions on the Occasion of His 70th Birthday. Editors: I. Hubeny, J. M. 
Stone, K. McGregor and K. Werner. American Institute of Physics. AIP Conference Proceedings 1171, 
pp. 27 -42, (2009). 



In contrast, the non-LTE problem of the 2nd kind is still rather unfamiliar to as- 
trophysicists. This term was introduced by Landi DeglTnnocenti |4|] to refer to the 
formidable numerical problem that implies calculating, at each spatial point of the (gen- 
erally magnetized) astrophysical plasma model under consideration, the values of the 
(27+ l) 2 elements of the atomic density matrix corresponding to each atomic level of 
total angular momentum J, which quantify its overall population, the population imbal- 
ances between its magnetic sublevels, and the quantum coherences between each pair 
of them. The values of such density-matrix elements have to be consistent with the in- 
tensity, polarization and symmetry properties of the radiation field generated within the 
medium. This requires solving jointly the RT equations for the Stokes parameters and 
the statistical equilibrium equations for the elements of the atomic density matrixH. Al- 
though much remains to be done, the numerical solution of the non-LTE problem of the 
2nd kind is facilitating the physical understanding and modeling of the so-called second 
solar spectrum, as well as the development of the field of astrophysical spectropolarime- 
try. 

As clarified below, the second solar spectrum is the observational signature of ra- 
diatively induced population imbalances and quantum coherences in the atoms and 
molecules of the solar atmosphere. Magnetic fields produce partial decoherence via the 
Hanle effect, giving rise to fascinating observable effects in the emergent spectral line 
polarization. Interestingly, these effects allow us to "see" magnetic fields to which the 
Zeeman effect is blind within the limitations of the available instrumentation [e.g.,^. 
In fact, the Zeeman effect is of limited practical interest for the exploration of magnetic 
fields in hot (e.g., chromospheric and coronal) plasmas because the polarization that it 
induces in a spectral line scales with the ratio between the Zeeman splitting and the 
Doppler width, and this ratio turns out to be particularly small in such plasmas. More- 
over, the polarization of the Zeeman effect as a diagnostic tool is blind to magnetic fields 
that are randomly oriented on scales too small to be resolved. We may thus define "the 
Sun's hidden magnetism" as all the magnetic fields of the extended solar atmosphere 
that are impossible to diagnose via the consideration of the Zeeman effect alone. 

2. THE SECOND STELLAR SPECTRUM 

Figure 1 contrasts the Fraunhofer spectrum (i.e., the visible intensity spectrum of the Sun 
with its multitude of absorption lines) with the so-called second solar spectrum, which is 
the linearly polarized spectrum observed by Stenflo and collaborators in "quiet" regions 
of the solar atmosphere [see [BL LzL [M Hi- The fractional linear polarization (Q/I) of the 
bottom panel shows a number of conspicuous peaks on top of an aparently noisy Q/I 



Note that the first Stokes parameter, /, is the specific intensity at a given wavelength, while Stokes 
Q represents the intensity difference between linear polarization parallel and perpendicular to a given 
reference direction (e.g., the direction perpendicular to the straight line joining the solar disk center 
with the observed point). Stokes U would be then the intensity difference between linear polarization 
at +45° and —45° with respect to the chosen reference direction, where positive angles are measured 
counterclockwise by an observer facing the radiation. Stokes V is the intensity difference between right- 
handed and left-handed circular polarization. 



5000 5500 8000 6500 

Wavelength [A] 



O .,.,,.!,, ,,.^, . 1 



5000 5500 6000 6500 

Wavelength [A] 

FIGURE 1. The visible Fraunhofer spectrum (displayed as the normalized intensity of the continuum) 
versus the second solar spectrum (displayed as the fractional linear polarization, Q/I). The positive 
direction for the Stokes Q parameter is the parallel to the observed solar limb. Extracted from [8]. 



background. Such Q/I peaks correspond to A6562.8 (H«), A5890 (Na I D 2 ), A5644.14 
(Ti I), A5565.48 (Ti i), A5535.51 (Ba i), A5167.32 (Mg i), A4758.12 (Ti I) and A4742.79 
(Ti I). For example, the left panels of Fig. 2 show an enlarged portion around the 
wavelengths of the sodium doublet. Similarly rich Q/I spectral structure is found in 
many other spectral regions, especially at wavelengths smaller than 5300 A (e.g., see 
the right panels of Fig. 2, which show the Q/I profiles produced by Y II, Cr I and Fe 
I). Clearly, the second solar spectrum has a structural richness that often exceeds that 
of the ordinary intensity spectrum on which most of our astrophysics is based. It is also 
of interest to note that the linear polarization amplitudes tend to increase towards the 
blue and near-UV spectral regions. For example, between 3910 A and 4630 A we find 
spectral lines with Q/I amplitudes of the order of 1%, such as Ca I A4227, Sr II A4078, 
Ba II A4554 and Sr I A4607 [see@]. 

As is well known, the polarization of the Zeeman effect is due to the wavelength 
shifts between the % (AM = M u -Mi = 0) and a^ r (AM = ±1) transition!! The typical 
observational signature of the circular polarization produced by the longitudinal Zeeman 
effect is an antisymmetric Stokes V(X) profile whose amplitude scales with the ratio, 
between the Zeeman splitting and the Doppler broadened line width. The linear 
polarization amplitudes of the transverse Zeeman effect scale instead as ^ 2 and its 




Such wavelength shifts are due to the presence of a magnetic field, which causes the atomic energy 
levels to split into different magnetic sublevels characterized by their magnetic quantum number M. 
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FIGURE 2. The two left panels are enlarged portions of Fig. 1 showing the complex Q/I structure 
across the sodium doublet, while the two right panels focus on a very interesting spectral region showing 
lines of Y II, Cr I and Fe I. 



typical observational signatures are symmetric Stokes Q(X) and U (A) profiles with their 
wing lobes of opposite sign to the line center one. 

As seen in the bottom panels of Fig. 2, the linear polarization signals of the second 
solar spectrum have nothing to do with the transverse Zeeman effect. This can also be 
seen in Fig. 3, which shows a high- sensitivity spectropolarimetric observation of the 
solar limb atmosphere in the strongest (8542 A) chromospheric line of the Ca II triplet. 
The observed Stokes V /I profiles are clearly caused by the longitudinal Zeeman effect, 
but the Stokes Q/I and U/I signals have Gaussian-like shapes at almost all spatial 
positions along the direction of the spectrograph's slit. Interestingly, while the Stokes 
Q/I signal changes its amplitude but remains always positive along that spatial direction, 
the sign of the Stokes U/I signal fluctuates. 

The shapes of the Q/I profiles of the second solar spectrum are not always Gaussian- 
like, however. A classification of the various Q/I shapes of atomic lines has been 
recently proposed by [11]. Strong resonance lines tend to show a triple peak structure, 
as is the case with the Na I D2 line shown in the bottom left panel of Figure 2. However, 
the shape of the observed Q/I and U/I profiles may also show spatial variability, as 
illustrated in Fig. 4 for the chromospheric K-line of Ca II. Note that in this near-UV line 
there is no measurable hint of the longitudinal Zeeman effect, in contrast with the case 
of the near-IR line of Fig. 3 whose line center Stokes Q/I and U/I signals also originate 
in the bulk of the solar chromosphere. 



3. THE PHYSICAL ORIGIN OF THE 2ND STELLAR SPECTRUM 



If the second solar spectrum is not caused by the Zeeman effect, then what could its 
physical origin be? 




FIGURE 3. Spectropolarimetric observation of the Ca II 8542 A line in a very quiet region close to 
the solar limb, obtained with the Zurich Imaging Polarimeter (ZIMPOL) at the Franco-Italian telescope 
THEMIS of the Observatorio del Teide (Tenerife; Spain). The reference direction for Stokes Q is the 
tangent to the closest limb. From 1 10]. 
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FIGURE 4. Spectropolarimetric observation of the K-line of Ca II in quiet regions of the solar atmo- 
sphere, obtained with ZIMPOL at the McMath Telescope in Kitt Peak Observatory (USA). Left panel: 
observation close to the solar limb. Right panel: disk center observation. From 11211 . 



3.1. Atomic level polarization and anisotropic radiation pumping 

In stellar atmospheres there is another fundamental mechanism, besides the Zeeman 
effect, producing polarization in spectral lines. There, where the emitted radiation can 
escape through the stellar surface, the atoms are illuminated by an anisotropic radiation 
field. The ensuing radiation pumping produces population imbalances among the mag- 
netic substates of the energy levels (that is, atomic level polarization), in such a way that 
the populations of substates with different values of \M\ are different. This is termed 
atomic level alignment. As a result, the emission process can generate linear polariza- 
tion in spectral lines without the need for a magnetic field. This is known as scattering 
line polarization [e.g.,E3|,EIl]- Moreover, radiation is also selectively absorbed when the 
lower level of the transition is polarized [15, ill 17, EH]. Thus, the medium becomes 



dichroic simply because the light itself has the chance of escaping from it. 



The fact that the illumination of the atoms in a stellar atmosphere is anisotropic 
is easy to understand if we consider the case of a plasma structure embedded in the 
optically thin outer layers (e.g., a solar coronal filament), because the incident radiation 
comes mainly from the underlying quiet photosphere and is contained within a cone 
of half aperture oc<90°, with the vertex centered on the point under consideration. The 
larger the height above the visible stellar "surface" the smaller a and the larger the 

anisotropy factor w = \/2Jq /Jq, where 7° = § ^1 ^ is the familiar mean intensity and 

Jq~§ ^§272^^ 2 ~^vQ' w * m f 1 = cos ^ being the cosine of the angle between the 
ray and the stellar radius vector (hereafter, the vertical direction). Neglecting the Q. 
dependence of the incident intensity / ^, and assuming that it is unpolarized, it is easy 
to find that w = [1 +cosa]cosa/2, which shows that in this case 0<w<l, with w = 1 
for the limiting case of a unidirectional unpolarized light beam that propagates along the 
vertical direction. 

The radiation field is also anisotropic within a stellar atmosphere itself (i.e., at heights 
where the overlying atmospheric plasma is not optically thin), but in this case w can be 
positive or negative. As shown in the right panel of Fig. 5, at such heights the radiation 
propagating outwards shows limb darkening (i.e., it is predominantly vertical) while 
that propagating inwards shows limb brightening (i.e., it is predominantly horizontal). 
Therefore, there is competition, because "vertical" rays (i.e., with > l/y/3) make 
positive contributions to w, while "horizontal" rays (i.e., with < l/y/3) make neg- 



ative contributions to w. As seen in Fig. 4 of [19], the larger the gradient of the source 
function the greater the anisotropy factor of the pumping radiation field, and the larger 
the amount of atomic level polarization. Therefore, at a given height in the inhomoge- 
neous solar photosphere, the anisotropy factor of the solar continuum radiation in the 
(granular) upflowing regions is significantly larger than in the (intergranular) downflow- 
ing plasma [see fig. 2 of l5|] . 

A useful expression to estimate the amplitude of the emergent fractional linear polar- 
ization is the following generalization of the Eddington-Barbier formula, which estab- 
lishes that the emergent Q/I at the line center of a sufficiently strong spectral line when 
observing along a line of sight specified by \x = cos0 is approximately given by [see 
Q: 

Q/I « ^=(l-^)[w^ Je ai(J u ) -w§ u o$(J,)], (1) 

(2) (2) 

where Wjjj and Wj j are numerical factors that depend on the angular momentum values 

(2) (2) 

(7) of the lower (/) and upper (u) levels of the line transition (e.g., w 10 = 1, w 0l = and 

wfi = —1/2), while <7q = pQ / pQ quantifies the fractional atomic alignment of the upper 
or lower level of the line transition under consideratior@- Note that the Oq values in Eq. 



4 For example, p# (J = 1 ) = (Ni + N + N- 1 ) / \/3 and pi {J = 1 ) = (Ni - 2N Q + N- 1 ) / \/6, where N\ , No 
and A^_i are the populations of the magnetic sublevels of a level with 7=1. 
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FIGURE 5. The emergent Q/I profiles (for a line of sight with ji =0.1) of three line transitions 
calculated in a model atmosphere with T = 6000 K and B=0G. All these Q/I signals are solely due to 
the atomic level polarization that results from the anisotropic illumination illustrated in the right panel. 
Left panel: assuming that the lower level is unpolarized. Middle panel: taking into account the full impact 
of lower-level polarization. The reference direction for Stokes Q is the parallel to the closest stellar limb. 
From lit 



(1) are those corresponding to the atmospheric height where the line center optical depth 
is unity along the line of sight. 

Consider the three line transitions of Fig. 5, and the corresponding emergent Q/I pro- 
files obtained by solving numerically the scattering polarization problem in an unmagne- 
tized model atmosphere assuming a two-level atomic model for each line independently. 
The left panel corresponds to calculations carried out assuming that the lower level is 
completely unpolarized, while the middle panel takes into account the full impact of 
lower-level polarization. Note that when the atomic polarization of the lower level is 
taken into account then the "null line" (i.e., the one which has unpolarized emission 
because it has J u = 0) shows the largest Q/I amplitude. Isn't it fascinating? I mean the 
fact that "zero-field" dichroism (i.e., differential absorption of polarization components 
in the absence of any significant Zeeman splitting) is a very efficient mechanism for 
producing linear polarization in the spectral lines that originate in a stellar atmosphere. 



3.2. The Hanle effect 

In order to understand what the Hanle effect is it is first necessary to clarify that in the 
general case where polarization phenomena are taken into account, the full description of 
an atomic system requires us to specify, for each 7- level, a density matrix with (27+ l) 2 



element^]. The diagonal ones, pj(M,M), quantify the populations of the individual 
sublevels and the non-diagonal ones, pj(M,M'), the quantum coherences between each 
pair of them. We say that the quantum coherence pj(M,M') is non-zero when the atomic 
wave function presents a well defined phase relationship between the pure quantum 
states \ JM) and \JM ). The law of transformation of the density-matrix under a rotation 
of the reference system chosen for the specification of its elements indicates that it is 
actually very common to find non-zero coherences when describing the excitation state 
of an atomic system under the influence of anisotropic radiative pumping [e.g., see [Till . 

Consider a reference system whose Z-axis (the quantization direction of total angu- 
lar momentum) is chosen along the direction of the applied magnetic field and /-levels 
whose sublevels are not affected by possible crossings and/or repulsions with the sub- 
levels pertaining to other levels. In this simplest case, the Hanle effect tends to reduce 
and dephase the quantum coherences with respect to the non-magnetic case, without 
modifying the population imbalances. For the Hanle effect to operate the magnetic field 
must be inclined with respect to the symmetry axis of the pumping radiation field. What 
happens is that as the sublevels are split by the magnetic field, the degeneracy of the 7- 
level under consideration is lifted and the quantum coherences are modified. This gives 
rise to a characteristic magnetic-field dependence of the linear polarization of the emer- 
gent spectral line radiation that provides an attractive diagnostic tool of magnetic fields 
in astrophysics [e.g., the review by[l9ll. 

Typically, in 90° scattering geometry (e.g., when observing off the solar limb) the 
largest polarization occurs for the unmagnetized case, with the direction of the linear po- 
larization perpendicular to the scattering plane. In the presence of a magnetic field point- 
ing either towards the observer (case a) or away from him/her (case b) the polarization 
amplitude is significantly reduced with respect to the unmagnetized case. Moreover, the 
direction of the linear polarization is rotated with respect to the zero field case. Normally, 
this rotation is counterclockwise for case (a) and clockwise for case (b). Therefore, when 
opposite magnetic polarities coexist within the spatio-temporal resolution element of the 
observation the direction of the linear polarization is like in the unmagnetized reference 
case, simply because the rotation effect cancels out. However, the polarization ampli- 
tude is indeed reduced with respect to the zero field reference case, which provides an 
"observable" that can be used to obtain information on hidden, tangled magnetic fields 
at subresolution scales in the solar atmosphere [see ll3l l5|l. 

On the other hand, in forward scattering geometry (e.g., when the line of sight points 
to solar disk center) we normally have zero polarization in the absence of magnetic 
fields, while the largest linear polarization (oriented parallel or perpendicular to the 
direction of the magnetic field) is found for a horizontal field with a strength such that 
the ensuing Zeeman splitting is sensibly larger than the level's natural width. In forward 
scattering geometry the linear polarization is created by the Hanle effect, a physical 
phenomenon that has been demonstrated to operate on the Sun via spectropolarimetric 
observations of coronal filaments in the He I 10830 A multiplet [seeQ/ZO. Another very 



5 If the atom under consideration has non-zero nuclear angular momentum then a rigorous modeling of the 
spectral line polarization requires taking into account the hyperfine structure, which drastically increases 
the number of density matrix elements required to specify the excitation state. 



interesting example is shown in the right panel of Figure 4. 

Approximately, the amplitude of the emergent spectral line polarization is sensitive 
to magnetic strengths between 0. 1 Bn and 10 Bu, where the critical Hanle field intensity 
(Bu, in gauss) is that for which the Zeeman splitting of the /-level under consideration 
is equal to its natural width: 

flH = 1.137xl(T7(f llfe gy) (2) 

with fiif e the lifetime, in seconds, of the J- level under consideration and gj its 
Lande factor. If the line's lower level is the ground level or a metastable level, its 
twe{Jl)~^/{BluJo)~l/{nAui)] (with n = c 2 J^/(2hv 3 ) being the solid-angle average 
of the number of photons per mode at the frequency of the line under consideration 
and B/ u and A u i the Einstein's coefficients). For relatively strong solar spectral lines 
?iif e (7/) is typically between a factor 10 2 and 10 3 larger than the upper-level lifetime 
(t]if & (J u )pdl / A u i, where A u \ is Einstein's coefficient for the spontaneous emission pro- 
cess). For this reason, in solar-like atmospheres the lower-level Hanle effect is normally 
sensitive to magnetic fields in the milligauss range, while the upper-level Hanle effect is 
sensitive to fields in the gauss range. 



4. THE NON-LTE PROBLEM OF THE 2ND KIND 

The generation and transfer of polarized radiation in a (generally magnetized) astro- 
physical plasma is a very involved non-local and non-linear RT problem which requires 
jointly solving the statistical equilibrium equations for the elements of the atomic den- 
sity matrix and the Stokes-vector transfer equation for each of the allowed transitions 
in the chosen multilevel model. Once the selfconsistent values of the elements of the 
atomic density matrix are known at each point within the medium, it is straightforward 
to solve the Stokes vector transfer equation in order to obtain the emergent Stokes pro- 
files. At present, there are two computer programs for solving non-LTE problems of 
the 2nd kind [see Both codes are based on the formal solvers and the iterative 



methods proposed by Trujillo Bueno for the numerical solution of non-LTE polarization 
transfer problems fsee 1 16L I3M. 



4.1. The Stokes-vector transfer equation 

In the polarized case, instead of the standard RT equation for the specific intensity 
7(v,n) one has to solve, in general, the following vectorial transfer equation for the 
Stokes vector I(v, 5) = (/,g,£7,V) t (f = transpose): 

^I=e-KI, (3) 

as 

where s measures the geometrical distance along the ray of direction Q., e is the emission 
vector, and K the propagation matrix. The elements of the emission vector (that is, 
£/, £q, £u, and £v) account for the contribution of the spontaneous emission process 



to the intensity and polarization that is generated at each spatial point by the physical 
mechanisms under consideration. The propagation matrix K can be written as 
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which helps to clarify that it has three contributions: absorption (the first matrix, Ki, 
which is responsible for the attenuation of the radiation beam irrespective of its polar- 
ization state), dichroism (the second matrix, K2, which accounts for a selective absorp- 
tion of the different polarization states), and dispersion (the third matrix, K3, which de- 
scribes the dephasing of the different polarization states as the radiation beam propagates 
through the medium). Actually, the situation is more complicated (and interesting!) be- 
cause in general the propagation matrix is the difference of two matrices: K = K A — K s , 
where K A is the above-mentioned contribution resulting from transitions from the lower- 
level towards the upper-level, while K s is the contribution caused by stimulated emission 
processes. Thus, Kj would be the amplification matrix (responsible for the amplification 
of the radiation beam irrespective of its polarization state), K| would be the dichroism 
matrix (responsible of a selective stimulated emission of different polarization states), 
and K| would be the corresponding dispersion matrix. Taking K| and K| into account 
may be of great interest for new discoveries in the field of astronomical masers [e.g., 

The general expressions of the components of the emission vector and of the propaga- 
tion matrix are very involved and will not be written here [see Chapter 7 of 14]. They are 
given in terms of the pQ multipolar components of the atomic density matrix (£ = 0,1,2 
and —K<Q<K) corresponding to the upper and lower levels of the line transition under 
consideration and of line-shape profiles, whose dependence on the magnetic quantum 
numbers cannot be neglected when the Zeeman splittings are a significant fraction of 
the spectral line width. Fortunately, such general expressions simplify considerably for 
several cases of practical interest . For instance, the following are the expressions of the 
only non-zero emission coefficients for the case of a one-dimensional atmosphere, either 
unmagnetized, or permeated by a microturbulent and isotropically distributed magnetic 
field of strength B, or in the presence of a magnetic field with a fixed orientation but with 
a strength in the saturation regime of the Hanle effect: 



£ e ne = £ o^-^=(i-M 2 )Po 2 , 



(5) 
(6) 



In this regime the quantum coherences turn out to vanish in a reference system whose quantization axis 
is parallel to the magnetic field vector (i.e. the real and imaginary parts of the and p| components are 
zero). 



where the pQ (K = 0, 2) values are those of the upper level of the line transition under 
consideration, £o = (hv /4n)A Ll £(j) x (2J u + X) x l 2 jV (with h, jY and (j) x being, respectively, 
the Planck constant, the particle density and the normalized line profile, with x the 

(2) 

frequency distance from the line center in units of the Doppler width), Wj j is the 
previously mentioned numerical factor that only depends on the quantum numbers of 
the levels involved in the transition, and where the orientation of the ray is specified 
by the azimuthal angle % and by \x = cosO (with 6 the angle between the direction 
of the radiation beam and the quantization axis, which must be the magnetic field 
direction itself for the deterministic field case and the normal to the surface of the stellar 
atmospheric model for the two other cases). Likewise, the only non-zero elements of 
the propagation matrix are r\[ and t\q which are given by identical expressions (i.e. by 

r\i = £/ and r\ Q = e Q ), but with r\ Q = (hv/4-7t)Bi u (j) x (2Ji + \) y l 2 jV instead of e , wfj u 

instead of xvf^j and with the p(f values of the lower level of the line transition instead 
of those of the upper level (for the case in which stimulated emissions are neglected). 
Note that £/ and Tjj depend on both the level's overall population (p°) and the population 
imbalances (Pq), while £q and r\Q depend only on the population imbalances (Pq A 

Under such assumptions a very useful formula can be derived to estimate the emergent 
fractional linear polarization amplitude at the center of a spectral line: 

Q_ g e /£ / -T7 g /T? / 

/~i-(i7eM)(£e/e/)' 

where £/, £q, rjj and t\q must be evaluated at height in the model atmosphere where 
the line-center optical depth is unity along the line of sight. This expression can be 
simplified taking into account the low polarization level in solar-like atmospheres (i.e., 
that r\Q/r\i -C 1 and £g/£/ 1). Thus, to first order in £g/£/ and r\Q/r\i, the fractional 
polarization reduces to <2//~£g/£/ — tiq/tji, which leads directly to the approximate 
Eq. (1). Its first term is due to selective emission of polarization components (caused 
by the population imbalances of the upper level), while its second term is due to "zero- 
field" dichroism (i.e., it is caused by the selective absorption of polarization components 
produced by the population imbalances of the lower level). Eq. (1) is very useful for 
understanding that, in general, the second solar spectrum is produced by such two 
contributions [see l23|l . 



4.2. The rate equations for the elements of the atomic density matrix 



We need also to consider the master equation for the atomic density matrix [see |24 , 
[HI]. Perhaps, the simplest and most famous example of this approach is that formulated 
by Einstein, who in 1917 introduced rate equations describing the effect of absorption, 
stimulated emission, and spontaneous emission processes in a two-level atom. The 



7 In these expressions, the reference direction for Q > is the perpendicular to the plane formed by the 
propagation direction and the Z-axis. 



situation is considerably more complex (and interesting!) in the present polarized case 
in which we have atomic and light polarization instead of simply overall populations of 
the atomic levels and intensity. In general, we have: 



^Pe(^) — [^PgC-OWagneticFields + t^Pg (^)i Radiation + [^Pg U)] Collisions' 

where the first two terms on the RHS are the contributions due to a magnetic field 
(e.g., via the Hanle effect) and to radiative transitions, while the last term is that due 
to collisions. In order to calculate the pg unknowns at each grid point of the chosen 

atmospheric model one usually assumes statistical equilibrium (i.e., dpg(/)/dt = 0). 

The quantum theory of spectral line polarization considered here assumes that the ab- 
sorption and re-emission processes are statistically independent events. When the only 
significant quantum coherences are those between the sublevels pertaining to each de- 
generate level, this complete redistribution theory [see [HI] provides a suitable descrip- 
tion if the pumping radiation has no spectral structure across a frequency interval larger 
than both the Larmor frequency (see below) and the natural width of the atomic levels. 

Of great interest is the Hanle-effect contribution of Eq. (8), which in the magnetic 
field reference frame reads [e.g., Eq. 7.11 in 1 1411 

[^Pe( 7 )]Hanle= " 2^ Yt*J GPfi GO, (9) 

where the Larmor frequency Vl = 1.3996xl0 6 5 (with B expressed in gauss) and gj 
is the Lande factor of the /-level under consideration. This expression, which is valid 
for the case of a multilevel atom assuming that quantum coherences can exist only 
between the magnetic sublevels pertaining to each individual /-level, indicates that in 
the magnetic field reference system the population imbalances (i.e. the pg density-matrix 
elements with 2 = 0) are insensitive to the intensity of the magnetic field, while the 
coherences (i.e. the Pg elements with Q^O) are reduced and dephased as the magnetic 
field strength increases. In fact, when the Zeeman splitting is much larger than the 
inverse lifetime of the shortest lived levels (e.g., when B > B satUT (J u ) rs 10Bh{J u )) such 
coherences vanish in the magnetic field reference system. 

Certainly isotropic collisions with electrons and ions contribute to the rate of change 
of the Pq(J) elements, but the most relevant effect in solar-like atmospheres comes 
typically from isotropic elastic collisions with neutral hydrogen atoms which induce 
transitions between the magnetic sublevels pertaining to the /-level under consideration. 
It can be shown that for isotropic elastic collisions one has [e.g., see Eq. 7.101 in 1 14f| 

[|pgWWcom S ions = -D^(J) P §(J), (10) 

where D^ K > (/) is the depolarizing rate (in s _1 ) whose value is proportional to the neutral 
hydrogen number density. Obviously A' ) (/) = 0, and the effect on the pg (/) elements 

with K > will be significant when 8^ (J) = D K (/) tu&(J) ss 1, where t U f e (J) is that 
defined after Eq. (2). 



Finally, the contribution caused by radiative transitions includes the transfer rates 
due to absorption (7^), spontaneous emission (7g) and stimulated emission (T$) from 
other levels, and the relaxation rates due to absorption (Ra), spontaneous emission (Re) 
and stimulated emission (R$) towards other levels. The general expressions are rather 
complicated [see, e.g., §7.2.aof[H, but the following ones are of great practical interest 
even though they neglect the stimulated emission rates and correspond to the particular 
case of a one-dimensional atmosphere, either unmagnetized or in the presence of a 
deterministic magnetic field with a fixed orientation but with a strength in the saturation 
regime of the Hanle effec1@: 

J, K, 

+ EE Po u ( J u) Wu,Ku -> J;K) 

~T,Po (J) [Ra(J;K,K -> /„) + R E (J;K,K -> /,)] , (11) 

k! 

where K x = 0, 2, 4, ... , 2J t ; K u = 0, 2, 4, 2J U and K* = 0, 2, 4, ... , 27. 

Of particular interest are the transfer rate 7a due to absorption from lower levels and 
the relaxation rate R& due to absorptions toward upper levels. Their expressions are given 
by linear combinations of the and radiation field tensors [see Eqs. 7. 14a and 7. 14d 
in [bill. Note that is the familar /-quantity of the standard non-LTE problem of the 
1st kind (which is an average over frequencies and directions of the Stokes / parameter 
weighted by the absorption profile), while! 

P = |dx^^^i=[(3 A i 2 -l)/^ + 3(m 2 -1)Q^]. (12) 

5. AN INTERESTING EFFECT: MAGNETIC MODULATION OF 

THE POPULATION RATIO 

Equations (11) imply that even the ratio jV u j 1 JV[ of the populations of the upper and 
lower levels of a spectral line may depend on the tensor component (which quantifies 
the "anisotropy" of the radiation field), and not only on Jq. Even more interesting is the 
fact that the jV u I jV\ value is modulated by the inclination of the magnetic field with 
respect to the symmetry axis of the incident radiation field (hereafter, the vertical axis). 
In order to illustrate this remarkable phenomenon let us consider the case of anisotropic 
radiation pumping in a gas composed by two-level atoms with J\ = 1 and J u = in 



It is important to note that in this last case the multipolar components, Pq (J), and the radiation field 
tensors, Jq, have to be defined in the magnetic field reference frame. 

9 In the most general case the number of Jq tensors necessary to quantify the symmetry properties of the 
radiation field is 9 [see §5.1 1 in 14] [see also Fig. 3 in l25tl . 



the presence of a deterministic magnetic field inclined by an angle 6b with respect to 
the vertical axis and with a strength B > 5 satur (7/)~10.B//(7/), so that all the quantum 
coherences between the magnetic sublevels of the lower level are zero in the magnetic 
field reference frame. In this Hanle effect saturation regime the statistical equilibrium 
equations in the magnetic field reference frame are identical to those corresponding to 
the unmagnetized case [see Eqs. 3-5 in ll5[], with the only difference that 7^ has to be 
substituted by: 

[%] B = ^(3cos 2 e fi -l)[^] v , (13) 

where [Jq]b an d [7q]v are the 7^ radiation field tensor in the magnetic field and in 
the "vertical" reference frames, respectively. Neglecting the contribution of inelastic 
collisions in Eq. (3) of 1150 one obtains 



where n = 7^(27/ + l)5/«/ (27 M + 1)A M / is the average of the number of photons per mode 
of the incident radiation at the line frequency and [w]b = v2[7q]b/7q is the anisotropy 
factor in the magnetic field reference frame. In Eq. (14) [Oq(7/)]b = [Pq{J)]b / Pq(J) is 
the fractional atomic alignment of the lower-level in the magnetic field reference frame, 
which can be easily calculated by neglecting inelastic collisions in Eq. (4) of [fl5ll . 
Disregarding also for simplicity elastic collisions and using the ensuing expression in 
Eq. (14) it is easy to find that for a two-level atom with 7/ = 1 and J u = 

^u = fir (3cos 2 9 B -l)V/2 1 

J/i 3 [ 4-(3cos 2 e B -l)w J ' 1 ' 



The lower curve in Fig. 6 shows the variation with Qg of jV u j 1 jY\ for this case of a 
resonance line transition with 7/ = 1 and J u = 0, assuming n = 10~ 4 (which implies that 
the rates of stimulated emission can indeed be neglected) and w = 0.99 (i.e., that the 
atomic system is irradiated by an almost collimated radiation beam£^- Note that when 
w = 1 the upper level is fully depopulated if the magnetic field is parallel to the incident 
beam of unpolarized radiation (i.e., if 6b = 0° or 6b = 180°). This occurs because, under 
such anisotropic pumping conditions, all atoms end up in the Mi = sublevel of the 
lower level and the medium becomes transparent. The same happens in the absence of 



magnetic fields, a well-known effect in the quantum optics literature [e.g.. 1261]. 

A similar but more involved analytical expression for JV u j ' jVi can also be derived 
for the case of a resonance line transition with 7/ = J u = 1, but assuming now that 
B > S satur (7 M )~10S//(7 M ). To this end it suffices with using the statistical equilibrium 
equations for the case of a two-level atom with 7/ = J u = 1 [seeded]. The upper curve 
in Fig. 6 shows the ensuing results. Note that for this case ^V u / = n = 10~ 4 when 



u It is of interest to point out that such values for n and w are very similar to those corresponding to the 
10830 A solar continuum radiation at a distance of 10 solar radii above the visible stellar surface. 
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FIGURE 6. Variation of the indicated population ratio with the inclination of the magnetic field vector 
with respect to the symmetry axis of the incident radiation. This radiation has h = 10~ 4 and w = [w]y = 
0.99. The strength of the magnetic field is assumed to be sufficiently high so as to produce saturation of 
the Hanle effect for the spectral line under consideration (e.g., only 0. 1 G for the line with J u =0 and 100 
G for the line with Ji=J u = 1). The dashed line indicates the n/3 value. 

[w]b = (which is always the case at B = 54.74° and 9 B = 125.26°, the so-called Van- 
Vleck angles). 

As seen in Fig. 6, for both line transitions the population ratio jV u / ' jY[ depends on 
the inclination of the magnetic field with respect to the symmetry axis of the pumping 
radiation field. Note that for the operation of this interesting effect in astrophysical 
environments where the stimulated emission terms are negligible (i.e., because n-Cl) 
the lower level of the line transition under consideration must be polarized. However, if 
stimulated emission radiation is significant (i.e., because n~ 1) then such a modulation 
takes place even in line transitions whose lower-level is unpolarized, such as in a spectral 
line with 7/ = and J u = 1 . 

6. CONCLUDING COMMENTS 

As mentioned in § 1 , we may define the Sun 's hidden magnetism as all the magnetic fields 
of the extended solar atmosphere that are impossible to diagnose via the consideration 
of the Zeeman effect alone. Contrary to what one might think, there are many examples 
that belong to this category: 

• The greater part of the magnetism of the quiet solar photosphere. 

• The magnetic fields of the solar chromosphere outside sunspots. 

• The magnetic fields that confine the plasma of solar prominences and filaments. 

• The magnetism of the solar transition region and corona. 

Fortunately, the Hanle and Zeeman effects can be suitably complemented for explor- 
ing magnetic fields in solar and stellar physics. To this end, it is necessary to interpret 



high-sensitivity spectropolarimetric observations through the numerical solution of non- 
LTE problems of the 2nd kind. The interested reader will find abundant information 
on a variety of interesting applications in solar physics in the following review papers: 

Finally, it is important to emphasize that one of the greatest challenges in astrophysics 
is the empirical investigation of the magnetic field vector in a variety of astrophysi- 
cal systems, such as the solar corona, circumstellar envelopes, accreting systems, etc. 
The physical mechanisms I have considered here (anisotropic radiative pumping, atomic 
level polarization and the Hanle and Zeeman effects) occur in many astrophysical envi- 
ronments, not only in the solar atmosphere. In particular, the linear polarization signals 
produced by the presence of atomic level polarization are sensitive to magnetic fields 
in a parameter domain that ranges from very weak fields (e.g., 1 mG) to a few hundred 
gauss. The observation and physical interpretation of these polarization effects via the 
solution of the non-LTE problem of the 2nd kind provides key information, impossible 
to obtain via spectroscopy or conventional spectropolarimetry. 
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